Dispersion of the high-energy phonon modes in Ndi.ssCeo.isCuO^ 
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The dispersion of the high-energy phonon modes in the electron doped high-temperature super- 
conductor Ndi.85Ceo.i5Cu04 has been studied by inelastic neutron scattering. The frequencies of 
phonon modes with Cu-0 bond-stretching character drop abruptly when going from the Brillouin 
zone center along the [100]-direction; this dispersion is qualitatively similar to observations in the 
hole-doped cuprates. We also find a softening of the bond-stretching modes along the [110]-direction 
but which is weaker and exhibits a sinusoidal dispersion. The phonon anomalies are discussed in 
comparison to hole-doped cuprate superconductors and other metallic perovskites. 



I. INTRODUCTION 

The question about the role of electron phonon cou- 
pling in high-temperature superconductivity is still mat- 
ter of controversy 0, IE B @ an d, most likely, will not 
be solved in the nearest future. Therefore, it appears 
interesting to continue to look for direct evidence for 
electron phonon coupling in these materials. In the con- 
ventional superconductors Q signatures of strong elec- 
tron phonon coupling can be found either in the form 
of dips in the phonon dispersion or by comparison with 
a non-superconducting reference material. In the case 
of the cuprates, however, there are no metallic non- 
superconducting reference materials; therefore, a safer 
way to discern electron-phonon coupling effects is to com- 
pare the phonon dispersion curves of cuprate supercon- 
ductors to those of their insulating parent compounds. 
The latter approach is adopted in this study. 

Very early, evidence for strong electron phonon cou- 
pling in the cuprate high-temperature superconductors 
was found in inelastic neutron scattering (INS) exper- 
iments on La2- x Sr x Cu04 and YBa2Cu306+5 The 
comparison of the phonon density of states (PDOS) mea- 
sured on the insulating parent compounds, La2Cu04 and 
YBa2Cu306, with that obtained on the optimally 
doped superconducting materials, Lai.85Sro.i5Cu04 and 
YBa2Cu307, reveals a pronounced shift of spectral 
weight. Doping induces a frequency renormalisation of 
the modes with the highest energies which possess a lon- 
gitudinal bond-stretching character, see Fig. 1. In con- 
trast, the low and medium energy ranges in the PDOS 
exhibit less sensitivity to doping. More recently, similar 
effects were also reported for the electron doped super- 
conducting cuprates 0, E| . 

Studies of the phonon dispersion using single crystals 
and triple axis neutron spectrometry have clarified the 
shift of weight in the PDOS 10]. Along the [100] di- 
rection, the branches of the longitudinal bond-stretching 
character exhibit a softening when going from the zone- 



centre into the zone, for the polarization patterns of the 
bond-stretching modes see Fig. 1. Qualitatively, the 
dispersion is similar in Lai.85Sro.i5Cu04 0, 0| and 
in YBa 2 Cu30 7 E Q 0. The branches are nearly 
flat close to the Brillouin zone centre, then show a 
steep decrease in frequency around q=(0.25 0) and 
become flat again between (0.35 0) and (0.5 0). 
Along the [110] direction a nearly flat dispersion is ob- 
served in La2- x Sr :c Cu04 and in YBa2Cu307, which 
still represents sizeable frequency renormalisation com- 
pared to the insulating parent compounds. The strong 
frequency renormalization of the zone-boundary modes 
is, however, not a particularity of the superconduct- 
ing cuprates, but is seen in many doped perovskites : 
nickelates [lfl^3, manganates [l8j | and superconducting 
Bai _^K^BiQc. ]l^. Eol El| . the only metallic perovskite 
compound which does not exhibit the effect is supercon- 
ducting Sr 2 Ru0 4 |22. 

We note that that the softening of the longitudinal 
bond-stretching phonons cannot be understood from phe- 
nomenological lattice dynamical models like the shell 
model. The insulator-to-metal transition upon doping 
entails a certain frequency reduction as a result of screen- 
ing effects suppressing the LO-TO splitting of the polar 
modes at the Brillouin zone centre. However, in mate- 
rials like the cuprates screening effects are confined to 
the long-wavelength phonons as it is illustrated in Fig. 
2. Such behavior is indeed observed in the dispersion of 
the low-energy phonon branches with polar character in 
many oxide perovskites [2, 0, [HI EH • 

Until recently the analysis of the bond-stretching 
phonon dispersion in the superconducting cuprates was 
restricted to La2- x Sr ;E Cu04 and YBa2Cu306+5, since 
only for these compounds it has been possible to grow 
single crystals of sufficient size for inelastic neutron scat- 
tering However, with the recent advance in synchrotron 
sources, inelastic x-ray scattering became more com- 
petitive and was applied to Nd2- x Ce :c Cu04 [2J, J29], 
HgBaoCuO/i [2f^ and very recently to La2- x Sr^Cu04 |27j j. 
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FIG. 1: Polarization patterns of the Cu-0 bond-stretching 
phonon modes for different propagation vectors : the zone 
center mode (a), the longitudinal mode for q= (0.25,0,0) 
(b), the linear or half-breathing mode (c), the longitudinal 
mode for q= (0.33,0.33,0) (d), the planar breathing mode 
(e), and the quadrupolar mode of transverse character at 
q=(0.5,0.5,0). 




FIG. 3: Polarization patterns of the zone-center modes rel- 
evant for the Ai high-energy dispersion. Large blue balls 
correspond to oxygen atom, small red balls to Cu atoms and 
the large yellow ones to Nd/Ce-sites. 



LO (unscreened) 




q - wavevector 

FIG. 2: Illustration of the screening of polar modes in a metal 
with low carrier concentration, when the electronic suscepti- 
bility is described in a free electron model. 



plain the phenomenological model used to analyze the 
data. Then, we present results for the [100] direction 
and the [110] direction and discuss the polarization pat- 
terns of the different phonon branches. Finally, we com- 
pare the phonon anomalies found in Nd2- x Ce :c Cu04 to 
those observed in hole-doped cuprates and also to those 
observed in other metallic oxide perovskite materials. 



For HgBa2Cu04 Uchiyama et al. [26| find a disper- 
sion similar to that observed in YBa2Cu307 HQ, i.e. 
the same step-like downwards dispersion in the bond- 
stretching branch along [100]. In the electron-doped 
material Nd2- x Ce ;E Cu04 the high-energy dispersion is 
more complex since there is another longitudinal branch 
of the same symmetry which interacts with the bond- 
stretching branch, see below. d'Astuto et al. 0? EH 
report that the frequencies of the modes with Cu-0 bond- 
stretching character drop near q~(0.2 0) falling below 
those of an other branch. However, the intensity related 
to bond-stretching character could not safely be followed 
beyond (0.15 0) in the inelastic x-ray scattering exper- 
iment. This motivated us to investigate the high-energy 
phonon dispersion in Ndi.gsCeo.isCuC^ by inelastic neu- 
tron scattering after large superconducting single crystals 
of this material became finally available. 

The paper is organized as follows: the next Section de- 
scribes our experimental technique. Thereafter, we ex- 



II. EXPERIMENTAL 

Two single crystals of Ndi.85Ceo.i5Cu04 were grown 
by a floating zone technique using a mirror furnace at Ky- 
oto University. The crystals exhibit the superconducting 
transition at 25 K and samples obtained with this method 
were already used in many studies [3, EH • Inelastic neu- 
tron scattering experiments were performed on the IT 
spectrometer at the Orphee reactor using Cu(lll) and 
Cu(220) monochromators and a pyrolithic graphite (002) 
analyzer. Higher order contaminations were suppressed 
with a pyrolithic graphite filter in front of the analyzer. 
Double focusing was applied on the monochromator and 
on the analyzer sides. Two single crystals, each of about 
500mm 3 volume, were co-aligned in the [100]/ [010] scat- 
tering geometry and cooled in a closed He-cycle cryostate. 
An additional experiment was performed with one crystal 
in the [100]/ [001] orientation. 
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III. RESULTS AND DISCUSSION 

The polarization patterns of the Cu-0 plane-polarized 
bond-stretching modes are depicted in Fig. 1. For 
the longitudinal modes at the zone-boundaries there is 
an alternation of Cu-sites with short and long metal- 
oxygen distances and therefore, the modes at (0.5,0,0) 
and (0,5,0,5,0) are called half-breathing and breathing 
modes, respectively. These modes may strongly couple to 
the charges, since smaller metal-oxygen distances point 
to a higher metal valence 30]. Fig. 1 depicts also the 
polarization patterns of some other plane-polarized Cu- 
01 bond-stretching modes whith a regard to discussion 
which modes can be expected to couple strongly to charge 
fluctuations. To this end, loci of dynamical charge accu- 
mulation favored by the atomic displacements are high- 
lighted. 

III.l Phenomenological model 

The above mentioned anti-crossing phenomena - which 
are absent in the case of La2- a; Sr :E Cu04 because of a 
different structure - makes the use of a lattice dynam- 
ics model indispensable to correctly assign the various 
high-energy phonon modes in Nd2-^Ce :E Cu04. This in- 
teraction of the longitudinal optical branches has been 
seen in the x-ray studies 0, |23| and is predicted by cal- 
culations with the phenomenological model. The com- 
plication arises from the shift of the second oxygen 
site (01 labels the in-plane oxygen and 02 the out- 
of-plane one) from the apical position at (0 0.17) in 
Lao_ ;E Sr ;E Cu04 to (0.25 0.25 0.25) in Nd 2 - x Ce x C\i04 
|31j . where it is stronger bonded with the RE- ions. In 
consequence, the vibrations of 02 parallel to the Cu02- 
planes exhibit higher energies and may interact with the 
Cu-0 bond-stretching modes in Nd2-^Ce :c Cu04. In or- 
der to put the assignment of the various high-energy 
modes onto firm ground, we measured the dispersion of 
all high-energy branches of a particular symmetry and 
compared the measured scattering intensities with the 
model predictions. Since it was beyond the scope of our 
experiments to determine a full set of dispersion curves 
we did not develop a new model but adapted the one 
fitted to the phonon dispersion of the parent compound 
Nd2Cu04 which is described in reference |22[. Only a few 
changes were necessary in order to achieve a satisfying 
agreement with the experimentally obtained dispersion. 
The model parameters are given in Table IL the crystal 
structure data were taken from reference |2]|. We use a 
shell model with Born-Mayer potentials describing 
the repulsive interactions. The shell charges of Nd and 
Cu and the amplitude of the 0-0 Born-Mayer potential 
were slightly adapted. Further, we included a longitu- 
dinal force constant between Cu and the 02-oxygen as 
was done for the non-doped compound 32] to arrive at a 



ionic part potentials 
ion Z Y K pair A(eV) r (A) C (eV/A b 
Cu 1.64 1.4 2.0 Nd-0 2000 0.319 
Nd 2.30 -4.9 8.0 Cu-0 3950 0.228 

01 -1.56 -3.0 2.0 0-0 2650 0.284 -100 

02 -1.56 -3.0 2.0 

force constants (dyn/cm) 
pair F 

Cu-02 (long.) 13372 

lin.-breathing 93800 
plan. -breathing 69300 

TABLE I: Lattice dynamics model parameters for the de- 
scription of the phonon dispersion in Ndi.ssCeo.isCuO^ for 
the explanation of the parameters see text; Z,Y are given in 
electron charges; K in 10 6 dyn/cm. 

fully satisfactory description of the experimental disper- 
sion curves. 

The main lattice dynamical difference between the par- 
ent and the doped compounds concerns the screening of 
the Coulomb-interaction through the free charge carriers. 
We describe the electronic susceptibility by the Lindhard- 
function for a free electron gas taking exchange correc- 
tions into account. Such model invariably produces com- 
plete screening at long distances, which means that the 
splitting between longitudinal optic (LO) and transverse 
optic (TO) phonon frequencies will disappear at the zone 
center. However, when passing into the zone, the screen- 
ing has to occur at smaller distances where it will be less 
effective and the LO-TO splitting may partially recover 
yielding an increasing dispersion. This scenario describes 
perfectly the dispersion of phonon branches with polar 
character at intermediate energies in many doped per- 
ovskites [HI [HI El EH • In particular, it applies to the 
modes with the strongest polar character which are the 
bond-bending modes related to ferroelectricity. 

In the case of Nd2-^Ce :E Cu04 we obtain a screen- 
ing vector of k s = 0.42 A -1 in good agreement with the 
analysis of the inelastic x-ray scattering results 24] . The 
screening-vector is furthermore comparable to the one 
found in La 2 - a ,Sr :c Cu04, k s = 0.39A- 1 for x-0.1 H3- 
In view of the higher doping in our sample the screening 
in the electron doped system appears to be somewhat 
less effective than in the hole-doped compounds. How- 
ever, in our analysis the screening is mainly determined 
by the slope of the 02-branch. Therefore, we cannot ex- 
clude that the screening is somewhat more effective in 
the Cu-0 plane. 

As has been pointed out in the introduction, screening 
treated within the free electron model is unable to explain 
the downward dispersion of the longitudinal branches 
along the [100] and the [110] directions. In order to 
reproduce such a behavior - which is sometimes called 
over-screening - we included additional force constants 
lowering the frequencies of the Cu-Ol bond- stretching 
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modes at the zone boundary in the [100] (linear breathing 
term) or in the [110] direction (planar breathing term), 
respectively. The linear and planar breathing constants 
( see Table I) correspond to force constants between a 
particular Ol-site and the other Ol's connected with the 
same Cu-atom. A displacement of the first 01 towards 
the Cu-site results in forces on the other Ol's also di- 
rected towards the Cu-site. With the linear and planar 
breathing constant, the forces act on just the opposite 
and on all three surrounding Ol's, respectively. With 
these parameters a quantitative description of the dis- 
persion is achieved |24|. We note that such terms have 
been used previously to describe the phonon anomalies 
in hole-doped compounds (see, e.g., jig). 

The softening of the zone-boundary modes observed in 
the superconducting cuprates can be considered as being 
due to an over-screening of the Coulomb-potentials re- 
quiring a special electronic susceptibility j^HI^. Assum- 
ing that charges may fluctuate between different cation- 
sites on the frequency scale of the phonon modes one may 
explain the softening of the breathing modes 0, HE l36j ]. 
Ab initio calculations of the phonon dispersion for the hy- 
pothetic material CaCu02 |23 an d for YBa2Cus07 
agree with the soft half-breathing mode frequency rather 
well, but they also obtain a soft breathing mode at 
q=(0.5 0.5 0) and they do not find the step-like disper- 
sion j^S E3 • Although the breathing force constants in 
our model do not posses a clear physical meaning they 
appear to mimic the particular electronic susceptibility 
in the metallic cuprates jsij- 

Table II resumes the group theoretical analysis. The 
zone center modes can be separated according to the irre- 
ducible representations, whose the polarization patterns 
are given in Table II. Along the [100] (A) and along the 
[110] direction (£), modes split into four different repre- 
sentations, all compatibility relations are given in Table 
II. 



III. 2 Phonon dispersion in [100] direction 

In the [100] direction the branches corresponding to 
the Ai-representation start from the E u -, A g - and B g - 
modes at the zone center. The polarization patterns of 
the four modes with the highest energies are shown in 
Fig. 3. The patterns were calculated with the lattice 
dynamical model fitted to the measured dispersion, see 
above. In addition to the CuOl bond-stretching mode 
(E u ) at the highest energy, there is the CuOl bond- 
bending mode (E u ) at about lOTHz, the i^-mode at 
10.3 THz and the 02-mode (E u ) at 9.4 THz. The B g - 
mode corresponds to a motion of 02 parallel to the c- 
axis; its frequency is observed by Raman-scattering at 
10.0-10.3 THz [II El [111. ^n-modes are Infra-Red ac- 
tive; however, the reported measurements are not conclu- 
sive, possibly due to insufficient sample quality [111. E3|. 
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TABLE II: Upper part : Polarization schemes according to the 
crystal structure of Nd2-xCe x Cu04 for all T- modes and for 
the representations along the [100] and [110] directions (la- 
beled A and E respectively); the first line gives the positions 
of 7 atoms forming a primitive unit, the following lines show 
the displacements of these atoms. A letter at the i-position, 
signifies that this atom is moving along the i-direction, a sec- 
ond appearance of the same letter signifies that the second 
atom moves with the same amplitude ("-" denotes a phase 
shift) in the corresponding direction. Lower part : compati- 
bility relations along the [100] and [110] directions 



Fig. 3 shows also the displacement pattern of the highest 
Eg mode. This mode is the starting point of a branch 
with A4 symmetry which lies in the same frequency range 
as those associated with the Cu-Ol respectively Nd-02 
bond-stretching vibrations. Measurements aiming at the 
dispersion of the latter branches were performed in Bril- 
louin zones with zero structure factor for the A4-modes. 

Energy scans at Q-values with favorable structure fac- 
tors for Cu-Ol-bond- stretching phonons are shown in 
Fig. 21 Results are shown both for medium and high 
resolution set-ups. At the zone center, (3 10), the bond- 
stretching mode can be easily detected. For intermediate 
q-values, however, it is obscured by the strong signal of 
the Nd-02 vibrations, whose frequency rapidly increases 
with q whereas frequencies associated with the Cu-Ol- 
bond-stretching vibrations abruptly decrease. The steep 
frequency increase in the 02-branch and the frequency 
decrease of the bond-stretching modes render the mea- 
surement very sensitive to the resolution conditions. At 
Q=(3.15 0) the increasing branch is focused and at 
Q=(2.85 ) the decreasing one, see Fig. 0] The interac- 
tion of the 01 and the 02-branch occurs around q=(0.2 
0). At q=(0.3 0) the energy of the Cu-Ol-bond- 
stretching mode is already more than 2 THz below that 
of the 02-mode. In addition, the Cu-Ol-bond-stretching 
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FIG. 4: Scans across the longitudinal bond-stretching modes 
with A i -symmetry: a)-d) are obtained with high resolution 
conditions (Cu(220) monochromator) , e)-h) with medium res- 
olution (Cu(lll)-monochromator) and i) with low resolution 
(Cu(lll)-monochromator with a final energy of 7.37 THz) . 
Lines correspond to fits with Gaussians. 
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FIG. 5: Dispersion of the high-energy phonon branches in 
Ndi.85Ceo.i5Cu04 in [100] direction. Symbols denote the 
experimental data and lines frequencies calculated with the 
phenomenological model. 
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FIG. 6: Resolution corrected linewidths observed for plane- 
polarized bond-stretching modes in Ndi.ssCeo.isCuCU in the 
[100] direction at T = 10 K. 



mode is broadened. Therefore, it has a stronger weight 
in the lower resolution configuration. At q=(0.5 0) the 
bond-stretching mode is clearly observed at 12.5 THz. 
The dispersion of the Ai-modes is resumed in Fig. 

We find a strong broadening of the Cu-Ol-bond- 
stretching modes for q-values in the Brillouin-zone. The 
full widths at half maximum were obtained by fitting 
Gaussian profiles to the scans and by correcting for the 
resolution, see Fig. Note that the experimental res- 
olution for measurements with the Cu(220) monochro- 
mator is much smaller (0.6 THz) than the width of the 
observed phonon peak. In spite of the difficulties to sep- 
arate Cu-0 bond-stretching modes from other ones our 
data unambiguously shows that these modes are signif- 
icantly broadened, yielding another evidence for strong 
electron phonon coupling in Nd2-^Ce iC Cu04. 

Due to the body centered stacking of the Cu02-planes 
in Nd2-a?Ce :E Cu04, the reciprocal lattice vector Z=(l 
0) is not a Brillouin-zone centre but a zone bound- 
ary: the corresponding modes are characterized by a 
phase shift of the displacements of atoms at r=(0 0) 
and at (0.5 0.5 0.5). Therefore, the displacement pat- 
terns of the modes at the zone center and of those at Z 
just differ by a phase shift in neighboring Cu02 planes. 
For the bond- stretching modes, one may expect simi- 
lar phonon frequencies at V and at Z due to the weak- 
ness of the inter-plane couplings involved. Following the 
Cu-Ol bond-stretching signal beyond (0.5 0) we in- 
deed find the expected frequency increase. However, 
between (0.6 0) and (1 0) there is an interaction 
with another branch. The Ai -branch starting at the B g - 
mode exhibits an increasing dispersion. Beyond q=(0.6 
0) these modes strongly mix with the Cu-Ol-bond- 
stretching ones. The high-energy phonon dispersion ob- 
tained by measurements in different Brillouin-zones, in- 
cluding zones with a large component, are resumed in 
Fig. Only a few measurements were made in the low 
energy region because we did not expect strong doping- 
induced changes here. The zone center modes in the 
range 8.5 to 10.5 THz were investigated in some detail be- 
cause there is a strong interaction between Ai branches of 
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FIG. 7: Polarization patterns of zone-center and Z-point 
modes with strong polarization along the c-direction. Large 
blue balls correspond to oxygen atom, small red balls to Cu 
atoms and the large yellow ones to Nd/Ce-sites. 



different polarization starting from these modes. These 
interactions virtually vanish, however, for q> 0.1 because 
the dispersion of the interacting modes is very different. 

We now turn to the c-polarized high-energy modes. 
Their polarization patterns are shown in Fig. [7| at 
the zone center, there is an A^-mode near 15 THz and 
the Bg-mode at 10.3 THz already mentioned. The 02- 
atoms exhibit a polarization along the c direction in 
both modes, but the modes differ with respect to the 
phase between neighboring 02-sites. In the A^-mode, 
all 02's move in phase, whereas they move out of phase 
in the I^-mode, see Fig. [71 At the Z-point, too, there 
are two c-polarized 02-modes corresponding to the in- 
phase and out-of-phase movement of the 02-ions at the 
same z-level. The B g mode is connected through the 
A i -branches with the in-phase movement at higher fre- 
quency; note that there are some interactions with other 
branches. The A^-mode connects with the in-phase Z- 
mode through the A4-branches. Typical scans aiming at 
the frequencies of these phonon modes are shown in Fig. 
[HI Compared to the longitudinal plane-polarized modes 
these phonons are better denned, i.e. their line-widths 
are smaller. The in-phase Ai-mode at Z is of special 
interest since all 02-sites move in-phase, which means 




12 13 14 15 16 17 18 10 11 12 13 14 15 16 
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FIG. 8: Several scans across high-energy phonons with differ- 
ent characters as indicated in the labels. Symbols denote the 
raw intensities and lines correspond to fits with Gaussians. 



that all neighboring 02 's of a particular Cu02-layer come 
closer at the same time, see Fig. [7| . Therefore, this 02- 
vibration may induce strong charge fluctuations between 
the Cu02 planes; it is related in character with the half- 
breathing mode, and one might expect a strong coupling 
with an inter-layer charge modulation. The correspond- 
ing mode in La1.85Sro.15Cu.O4, labeled Ozz, indeed is 
rather anomalous with a large line width and a frequency 
renormalisation compared to the parent compound by 
several THz [l(j. This effect is interpreted by Falter et al. 
Q as due to strong coupling with the inter-layer plasmon 
mode. We find that the O^^-mode is not well defined in 
Ndi.85Ceo.i5Cu04 as well indicating a large intrinsic line 
width. Unfortunately, it is difficult to isolate the con- 
tributions of the c-polarized Z-point phonon from those 
with the in-plane polarization , but the observation of a 
three-peak structure obtained at several Z-points clearly 
proves that the c-polarized mode exhibits a high energy. 
There is strong evidence that the O^z-mode has a large 
line width. Compared to Lai.85Sr .i5CuO4 the Ozz~ 
mode in Nd1.85Ceo.15Cu.O4 appears to be less anoma- 
lous. 



III. 2 Phonon dispersion in [110] direction 

The phonon dispersion in the [110] direction, E, was 
studied, again with the focus on the bond-stretching 
branches. Typical scans are shown in Fig. [2l and the 
phonon dispersion derived from the data is drawn in Fig. 
ITU1 Again, the analysis is complicated by the interaction 
of the Cu-Ol with the Nd-02 bond-stretching branches. 
The scattering associated with the latter modes is in most 
configurations dominant. Nevertheless, we were able to 
unambigously determine the dispersion of the branches 
of interest. The following picture emerged: the branch 
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v (THz) v (THz) 



FIG. 9: Scans across the longitudinal a-d) and transverse e- 
h) bond-stretching modes with Si and S3 symmetry, respec- 
tively. All scans were performed using a Cu(lll) monochro- 
mator with a final energy of Ef=3.55 THz, except that in g) 
and h) where E f =7.S7 THz. . 
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FIG. 10: Dispersion of the high-energy phonon branches in 
Ndi.85Ceo.i5Cu04 in [110] direction. Symbols denote mea- 
sured frequencies and lines results of the lattics dynamical 
calculations (the broken line in the middle corresponds to the 
S 2 representation. 

quadrupolar mode can be understood from the fact that 
this vibration involves weaker force constants compared 
to the zone center bond-stretching mode, as the Ol-Ol 
distances vary only little. 



associated with Nd-02 vibrations disperses steeply up- 
wards from its starting point at v = 9.4 THz and becomes 
the highest one in frequency for q larger than (0.15 0.15 
0). At the zone boundary it reaches a frequency of 17.4 
THz, the highest in the entire dispersion. Raman exper- 
iments have observed a peak at just this frequency 
suggesting that it should be attributed to the (0.5 0.5 
0) LO 02-mode. This mode, which is not Raman ac- 
tive, may become so through the disorder induced by the 
Ce-doping. However, it remains astonishing that this dis- 
order induced Raman-peak is that strong. In contrast to 
the Nd-02 branch, the Cu-Ol-bond-stretching branch is 
found to exhibit a downward dispersion similar to the ob- 
servation along the [100] direction. However, the down- 
ward dispersion is less pronounced along [110] and has a 
more conventional shape, i.e. it is sine-like. 

The transverse Cu-Ol-bond-stretching branch of £3- 
symmetry also exhibits a steep frequency decrease when 
passing into the Brillouin-zone, see the right part of Fig. 
EH and Fig. The transverse bond-stretching branch 
along the [110] direction ends at the zone-boundary in 
the so-called quadrupolar mode, whose polarization pat- 
tern is given in Fig. 1. Its low frequency compared to 
that of the zone centre mode is perfectly described by the 
lattice dynamical model without any need for a special 
force constant. The shape of this branch has to be con- 
sidered as normal. A similar shape is observed in many 
other perovskites. The relatively low frequency of the 



IV. Discussion of bond-stretching phonon anomalies 

- Comparision with inelastic x-ray scattering data - 
Qualitatively, our results agree with the anticrossing sce- 
nario of Cu-Ol bond-stretching and Nd-02 modes which 
was deduced from an inelastic x-ray experiment 0, [2^ , 
but quantitatively, there are significant differences for q- 
values away from the Brilluoin-zone center. For the [100] 
direction we find good agreement up to q=(0.15 0), 
whereas already at q=(0.2 0) the tentative identifica- 
tion of the bond- stretching mode near 12 THz seems to 
be incorrect, see Fig. [H] • The inelastic x-ray spectra 
are dominated by the contribution of the Nd-02 peak, 
therefore, the bond-stretching frequencies could not be 
determined unambiguously but had to be fitted to the 
little structured back ground. For instance, inspection of 
Fig. 2 in 24] shows that no clear peak was observed at 
the position of the Cu-Ol mode for q=(0.4 0). The ap- 
proximate agreement between the x-ray and the neutron 
frequencies for this q appears to be somewhat accidental. 

The problems encountered in the x-ray experiment 
aiming at the dispersion of the Cu-Ol bond-stretching 
modes in the [110] direction are similar to those just dis- 
cussed for the [100] direction. In particular, the inelastic 
structure factors of these modes are very small for q be- 
yond (0.25,0.25,0) in the Brillouin zone chosen for the 
measurements. Therefore, it is not surprising that the 
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x-ray data for this q-region are quite inaccurate and fall 
1.5 THz below the neutron data. 

In summary, our neutron results confirm the anti- 
crossing scenario discussed in 0> HE) but only qualita- 
tively so. The neutron experiment is much better suited 
to resolve the complex multi-branch dispersion by exam- 
ining many different Brillouin zones and due to the fact 
that oxygen and heavy-ion modes posses similar dynami- 
cal structure factors for neutron scattering. Higher order 
scattering events may play a stronger role in the case of 
inelastic x-ray scattering studies on the oxygen modes. 

- Comparision of bond- stretching phonon anomalies in 
high temperature superconductors - Due to the mixing of 
several branches in Ndi.ssCeo.isCuC^ it is not obvious 
to isolate the behavior of the bond-stretching modes in 
the [100] direction from Fig. [5J Therefore, we present 
the Ai- data again in Fig. EI where the frequencies of 
the modes with bond-stretching character are drawn as 
filled and all other contributions as open circles, respec- 
tively. The separation can be easily done in the first half 
of the Brillouin-zone, but is less obvious in the second 
half due to the interaction with the c-polarized modes, 
which is more extended in q-space. The presentation 
in Fig. El clearly shows the resemblance of the bond- 
stretching dispersion in Ndi.ssCeo.isCuC^ with that in 
Lai.85Sro.i5Cu04 |lll |. The observed dispersion is rather 
well reproduced by the shell model after inclusion of the 
additional breathing constants except that the frequency 
drop around q= (0.25,0,0) is sharper in experiment . We 
emphasize that the classical shell model without such 
force-constants predicts an upwards bond-stretching dis- 
persion as is indicated in Fig. El by broken lines. In the 
top of Fig. El we show the differences between the mea- 
sured bond-stretching frequencies and those calculated 
from the shell model without breathing constants. These 
differences demonstrate the influence of the anomalous 
strong electron-phonon coupling in Nd2- £C Ce £C Cu04. 

The anomalous bond-stretching dispersion in 
Ndi.85Ceo.i5Cu04 is compared with that of other 
cuprate superconductors in Fig. ^1 Evidently, the dis- 
persion in the [100] direction is step-like in all optimally 
doped cuprate superconductors investigated so far. The 
different levels of bond-stretching phonon frequencies 
can be attributed to different Cu-Ol-bond distances 
induced either through different dopants or through 
internal strain. When comparing the frequencies of the 
half-breathing mode to those of the zone-centre bond- 
stretching mode we find the strongest renormalisation in 
the two high-T c materials YBa2Cus07 and HgBa2Cu04; 
Lai.85Sro.i5Cu04 exhibits a slightly stronger softening 
than Ndi.85Ceo.i5Cu04 along [100] |4^. In contrast, the 
downward dispersion along the [110] direction is most 
pronounced in Ndi.85 Ceo. 15 CuC>4 whereas the corre- 
sponding branch in Lai.ssSro.isCuC^ is nearly flat and 
exhibits even a weak frequency increase in YBa2Cus07 
fTo| . That is to say, Ndi.ssCeo.isCuC^ shows the 



least anisotropy between the frequency renormalisation 
along the [100] and the [110] directions of all optimally 
doped cuprates studied so far. In this sense, it can be 
said that the phonon softening in Ndi.ssCeo.isCuCU 
has the least one-dimensional character. It remains 
to be seen whether this fact is related to the distinct 
character of the charge carriers, electron versus hole-like, 
in the cuprate high-temperature superconductors. It is 
worth mentioning that the one-dimensional nature of 
the over-screening effect is most pronounced in the 90 
K-superconductor YBa2Cu306- But we want to em- 
phasize once more, that the anomalous bond-stretching 
phonon dispersion in all superconducting cuprates is 
astonishingly similar. 

There is an on-going debate about the importance of 
electron phonon coupling for the kink-like features ob- 
served in ARPES spectra of many cuprate high temper- 
ature superconductors 0,0,0- However, there seem to 
be significant differences between the ARPES features of 
hole and electron doped compounds |43 , whereas we find 
the phonon anomalies in the bond-stretching branches in 
all cuprates to be essentially the same. This raises doubts 
about a dominant role of the bond-stretching modes for 
the ARPES features. 

-Comparison with other metallic perovskites - The 
over-screening effects in the bond-stretching dispersion 
are not a particularity of the superconducting cuprates 
but are seen in many doped perovskites : nickelates [lfil 

tmanganates [HI and superconducting Bai-^K^BiOs 
I20L El| , the only metallic perovskite compound which 
does not exhibit the effect is superconducting Sr2Ru04 
allthough this material is the best metal. The fre- 
quency softening is even much more pronounced in man- 
ganates and bismuthates than in the cuprates |2lj . 

Fig. El further demonstrates that many perovskite 
materials exhibit even a similar shape of the bond- 
stretching dispersion along [100]. The step-like shape 
of the dispersion is observed in Bai-^K^BiOa and in 
Lai-^Sr^MnOa for x=0.2 and 0.3 . In the nickelates, 
however, the dispersion is nearly sinusoidal. This dis- 
persion may be related to the insulating electronic be- 
havior in nickelates 0, The absolute size of the 
over-screening effects differs essentially between the sys- 
tems. The effect along [100] is strongest in the man- 
ganates where it amounts to a reduction by 30% . Along 
the other directions in reciprocal space, the manganates 
and Bai-^K^BiOs exhibit over-screening effects which 
even exceed those occurring along [100]. For instance, 
the volume breathing mode in Bai-^K^BiOa is renor- 
malized by 40% compared to the insulating reference ma- 
terial [23,1^. A strong softening in the [110] direction 
was observed also in the nickelates. This indicates that 
the electron-lattice coupling is quite strong in these ma- 
terials, even considerably stronger than in the cuprates. 
This conclusion is corroborated by the fact that in the 
manganites and in Bai-^K^BiOs the phonon anomalies 
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FIG. 11: (Below) Dispersion of the bond-stretching phonon 
frequencies (filled circles) and other modes (open circles) com- 
pared to the calculated phonon frequencies. The dashed and 
continous lines correspond to screened shell-model calcula- 
tions without and with the additional breathing-force con- 
stants, respectively. (Above) Difference between the mea- 
sured bond-stretching frequencies and the shell-model calcu- 
lations without the breathing forces. 

are not restricted to a line in reciprocal space but are very 
strong also in other directions. Therefore, they entail a 
much stronger shift of spectral weight in the PDOS. 

We think that there is a close relationship of the bond- 
stretching phonon anomalies with the charge ordering 
phenomena in each particular system [2l| . All these sys- 
tems including the cuprates are close to charge ordering 
instabilities, which may be stabilized by a particular dop- 
ing |43,|43,|43. Roughly speaking, those bond-stretching 
modes are most renormalized in frequency whose polar- 
ization patterns corresponds to the structural distortion 
induced by the charge ordering. The most obvious dif- 
ference between the cuprates and other perovskite com- 
pounds consists in the fact that the [100] direction plays 
a special role only in the cuprates, what perfectly reflects 
the one-dimensional nature of charge ordering in the form 
of stripes in hole-doped cuprates. In contrast, charge 
ordering phenomena in nickelates, manganates and bis- 
muthates occur with a checker-board arrangement. Con- 
currently, the strongest bond-stretching phonon anoma- 
lies in these compounds occur along the related [110] 
or [lll]-directions. The relation between charge order- 
ing and bond-stretching phonon anomalies also holds for 
the ruthenate, since in layered ruthenates, which exhibit 
a normal bond-stretching phonon dispersion, in-plane 
charge ordering seems not to occur. 

Charge ordering phenomena have so far not been ob- 
served in the electron-doped cuprates. From their similar 
phonon dispersion, one may expect them to appear most 
likely again along the [100]-direction. 

Anomalous peak splitting has been reported for a num- 
ber of perovskite materials and have been taken as fin- 
gerprints of an inhomogeneous charge distribution, i.e. in 




FIG. 12: Comparison of the phonon anomaly in the bond- 
stretching branches observed in different metallic oxide per- 
ovskite materials. In the left we show the dispersion along 
the [100] direction and in the right part that for the [110] di- 
rection. Data were taken from references : Unstoichiometric 
La 2 Ni0 4 |3 JLai.69Sro.3iNi0 4 17], Lai.ssSro.isCuCU E3, 
YBa 2 Cu 3 6 0^HgBa 2 CuO4 EfTBai-^BiOs EH, and 
Lai-.Sr^MnOs El 



Bai-^K^BiOs m the nickelate s Ull . and also in 

YBa 2 Cu 3 6 and Lai. 85 Sro.i5Cu0 4 0,0. For the lat- 
ter two materials, however, it was shown that there are 
no indications for intrinsic splittings of bond-stretching 
branches 0. We note that in Ndi.ssCeo.isCuC^ 
as well, all observed peak splittings of bond-stretching 
phonons can be fully explained with standard theory of 
lattice dynamics, i.e. due to the interaction of branches 
having the same symmetry. 



IV. CONCLUSION 

Inelastic neutron scattering allowed us to quantita- 
tively characterize the high-energy phonon dispersion in 
Ndi.85Ceo.i5Cu04 . Previous inelastic x-ray measure- 
ments had already shown that the frequencies of the 
Cu-0 bond-stretching branches are strongly renormal- 
ized when compared to those of the insulating parent 
compounds but reliable data could be obtained only for 
a limited range of wave vectors very close to the zone 
center whereas the frequencies of the linear and of the 
planar breathing modes were not accessible. 

The dispersion of the branches with Cu-Ol bond- 
stretching character exhibits anomalies astonishingly 
similar to the observations in hole doped cuprate su- 
perconductors. Along the [100] direction the frequency 
of the bond-stretching modes abruptly decreases near 
q=(0.25 0). This effect is slightly weaker than in the 
hole doped materials. Moreover, these phonons acquire 
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a large linewidth which underpins their anomalous char- 
acter. The longitudinal bond-stretching modes exhibit a 
downward dispersion also along the [110] direction with a 
frequency drop between zone center and zone boundary 
of about 1 THz but the dispersion shape of the branch 
appears to be more normal. The softening of the pla- 
nar breathing-mode softening in Ndi.ssCeo.isCuCU is 
stronger than the corresponding effect in the hole doped 
materials, in particular compared to YBa2Cu306, indi- 
cating that the anomalous softening of the linear breath- 
ing mode has weaker one-dimensional character. How- 
ever, the differences in the phonon anomalies between 
hole and electron doped cuprates concern details; the 
signatures of electron phonon coupling in the bond- 
stretching phonons are found to be essentially the same 
in all cuprate high temperature superconductors studied 
so far. 
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